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THE RATIONAL MOTION OF MINIMAL DUAL QUATERNION 
DEGREE WITH PRESCRIBED TRAJECTORY 


ZIJIA LI, JOSEF SCHICHO, AND HANS-PETER SCHROCKER 


Abstract. We give a constructive proof for the existence of a unique rational 
motion of minimal degree in the dual quaternion model of Euclidean displace¬ 
ments with a given rational parametric curve as trajectory. The minimal motion 
degree equals the trajectory’s degree minus its circularity. Hence, it is lower 
than the degree of a trivial curvilinear translation for circular curves. 


1. Introduction 

A rational motion is a motion with only rational trajectories. In the dual 
quaternion model of SE(3), the group of rigid body displacements, [12, Ch. 9] it is 
described by a rational curve on the Study quadric [6]. In this article we construct a 
rational motion of minimal degree in the dual quaternion model with a given rational 
curve as trajectory, and we show that this motion is unique up to coordinate changes. 
This is an interesting result in its own right but it also has a certain potential for 
applications in computer graphics, computer aided design or mechanism science. 

Usually, one defines the degree of a rational motion as the maximal degree of a 
trajectory [6]. With this concept of motion degree, our problem becomes trivial as 
the curvilinear translation along the curve is already minimal. As we shall see, it is 
also minimal with respect to the dual quaternion degree if the prescribed trajectory 
is generic. The situation changes, however, if the trajectory is circular, that is, 
it intersects the absolute circle at infinity. In this case, the minimal achievable 
degree in the dual quaternion model is the curve degree minus half the number 
of conjugate complex intersection points with the absolute circle at infinity (the 
curve’s circularity). 

We will see that twice the circularity of a trajectory equals the trajectory degree 
minus the degree defect in the spherical component of the minimal motion. This 
leads to the rather strange observation that generic rational motions (without 
spherical degree defect) have very special (entirely circular) trajectories. Conversely, 
the minimal motion to generic (non-circular) curves are curvilinear translations 
which are special in the sense that their spherical degree defect is maximal. 

We continue this article with an introduction to the dual quaternion model of 
rigid body displacements in Section 2. There we also introduce motion polynomials 
and their relation to rational motions. Our results are formulated and proved in 
Section 3. The proof of the central result (Theorem 2) is constructive and can be 
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used to actually compute the minimal rational motion by a variant of the Euclidean 
algorithm. We illustrate this procedure by two examples. 

2. The dual quaternion model of rigid body displacements 

In this article, we work in the dual quaternion model of the group of rigid body 
displacements. This model requires a minimal number of parameters while retaining 
a bilinear composition law. Moreover, it provides a rich algebraic and geometric 
structure. It is, for example, possible to use a variant of the Euclidean algorithm 
for computing the greatest common divisor (gcd) of two polynomials. This section 
presents the necessary theoretical background on dual quaternions. 

2.1. Dual quaternions. The set DH of dual quaternions is an eight-dimensional 
associative algebra over the real numbers. It is generated by the base elements 

1, i, j, k, £, d, £j, £k 

and the non-commutative multiplication is determined by the relations 

= k^ = ijk = —1, £^ = 0, i£ = £i, j£ = £j, k£ = £k. 

As important sub-algebras, the algebra of dual quaternions contains the real numbers 
K = (1), the complex numbers C = (l,i), the dual numbers D = (1,£), and the 
quaternions El = (1, i, j, k) (angled brackets denote a linear span). A dual quaternion 
may be written as h = p + eq where p, q are quaternions. The conjugate dual 
quaternion is h = p + eq and conjugation of quaternions is done by multiplying 
the coefficients of i, j, and k with —1. It can readily be verified that the dual 
quaternion norm, defined as ||/i|| = hh, equals pp + e{pq + qp). It is a dual number. 
The non-invertible dual quaternions h = p + eq are precisely those with vanishing 
primal part p = 0. 

An important application of dual quaternions is the modelling of rigid body dis¬ 
placements. The group of dual quaternions of unit norm modulo {±1} is isomorphic 
to SE(3), the group of rigid body displacements. A unit dual quaternion h = p + eq 
acts on a point x in the three dimensional real vector space (i,j,k) according to 

(1) X i-P- pxp + pq — qp = pxp + 2pq. 

Note that —qp = pq because of the unit norm condition. It is convenient and 
customary to projectivise the space DH of dual quaternions, thus arriving at P^, 
the real projective space of dimension seven. Then, the unit norm condition can 
be relaxed to the non-vanishing of pp and the vanishing of pq + qp. In a geometric 
language, SE(3) is isomorphic to the points of a quadric S C P^, defined by 
pq + qp = 0, minus the points of a three-dimensional space, defined by p = 0. The 
quadric S is called the Study quadric. In this setting, the map (I) becomes 

pxp + pq — qp pxp + 2pq 

PP PP ' 

The action oi h = p + eq with p ^ 0, pq + qp = 0 can be extended to real projective 
three-space P^, modelled as projective space over (1, i,j,k). The point x represented 
by xo + x\i + X 2 ^ + a^ak is mapped according to 

X !->■ XqPP + p{xii + X2j + X 3 k)p + 2xopq = pxp + 2xopq. 

This is a convenient representation for studying rational curves as trajectories of 
rational motions. 
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2.2. Rational motions and motion polynomials. In the projective setting, a 
rational motion is simply a curve in the Study quadric S that admits a parameteri- 
sation by a polynomial 

n 

(2) C = 

i=0 

with dual quaternion coefficients co,...,c„ € DM. The non-commutativity of 
DM necessitates some rules concerning notation and multiplication: Polynomial 
multiplication is defined by the convention that the indeterminate commutes with 
all coefficients, the ring of these polynomials in the indeterminate t is denoted by 
DM[t], the subring of polynomials with quaternion coefficients is M[t]. We always 
write coefficients to the left of the indeterminate t. This convention is sometimes 
captured in the term “left polynomial” but we will just speak of a “polynomial”. 

Evaluating C at different values t G M gives points of a rational curve in P^. We 
also define C(oo) := c„ = limt_>oo in order to obtain the complete curve as 

image of P^ := K U {oo} under the map 1 1 —> C'(t). A similar convention is used for 
rational curves in P^. 

The conjugate to the polynomial (2) is 

n 

C-.= Y,Cit\ 

i=0 

It can readily be verified that the norm polynomial CC has dual numbers as 
coefficients. If CC has even real coefficients and the leading coefficient c„ is 
invertible, we call C a motion polynomial. This is motivated by the observation that 
the curve parameterised by a motion polynomial is contained in the Study quadric, 
that is, it parameterises indeed a rigid body motion. Writing C = P + eQ with P, 
Q G M[t], motion polynomials are characterised by the vanishing of PQ + QP. 

The trajectory of the point x = xq + xj -I- X 2 i + xsk G P^ is 

(3) t !—>■ PxP + 2xoPQ, t G P^. 

This is a curve parameterised by polynomial functions in projective coordinates. It 
is also called a rational curve because it is always possible to clear denominators 
of rational functions. All motions in SE(3) with only rational trajectories can be 
parameterised by motion polynomials [6]. 

3. Rational curves as trajectories of rational motions 

A rational parameterised equation x = xq + xii -1- a: 2 j + with xq, Xi, X 2 , 
xs G M[t] is called reduced if the greatest common divisor g of xq, xi, X 2 , X 3 has 
degree zero. The degree of x is defined as the maximal degree of xq, xi, X 2 , X 3 . If x 
is not reduced, we may divide it by g to obtain an equivalent parameterised equation 
which describes the same rational curve, but possibly with fewer parameterisation 
singularities. 

Given a reduced parameterised equation x of degree d, it is our ultimate aim to 
find a motion polynomial C = P + eQ G DM[t] of minimal degree such that the 
trajectory of one point p is parameterised by x. An obvious example of a motion 
polynomial with trajectory x is the curvilinear translation along x and it will turn 
out that this is already the solution to our problem in generic cases. However, for 
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trajectories which are non-generic, in a sense to be made precise below, a lesser 
degree can be achieved. 

3.1. Circularity of trajectories. A motion polynomial C = P + eQ is called 
reduced, if both P and Q have no common real factor. The real factor of the primal 
part P with maximal degree is uniquely defined up to a constant scalar factor. We 
call it maximal real polynomial factor and abbreviate it by “mrpf”. It accounts for a 
difference in degrees between the rational motion, parameterized by C = P -I- eQ, 
and its spherical component, parameterised by P. Hence, we also call the degree of 
the maximal real polynomial factor of P the spherical degree defect of C. 

We only consider reduced rational motions but do not exclude motions with 
positive spherical degree defect. It will turn out that the spherical degree defect 
and circularity of trajectories are closely related: 

Definition 1. The circularity of a rational curve with reduced rational parame¬ 
terised equation x = xo-l-a;ii-l- X 2 j + is defined as 

^ deggcd(xo ,x\ + xl + xf). 

The curve is called entirely circular if it is of maximal circularity ^ degx. 

Geometrically, circularity is half the number of (conjugate complex) intersection 
points of X and the absolute circle at infinity, counted with their respective algebraic 
multiplicities. Hence, it is always an integer (this also follows algebraically from 
the fact that Xq, Xi, X 2 , X 3 are real and prime) and does not depend on the chosen 
parameterisation, as long as it is reduced. 

Theorem 1. Let C be a reduced motion polynomial of degree n and spherical 
degree defect m. Then a trajectory of C is of degree d < 2n — m and circularity 
c> {d — m)l 2 . 

Remark 1. In Theorem I, the strictly less and strictly greater cases can occur. One 
example of this is the general Darboux motion [7, 9] where n=3, m = d=2, c=0 
and hence 2 = d <2n — m = A (but 0 = c = {d — m)/2). Another example is the 
circular translation motion [7, 9] with n = m = d = 2, c = 1. Here 2 = d = 2n — m 
but 1 = c > (d — m)/2 = 0. 

Lemma 1. Assume A,Bg If a monic quadratic irreducible real polynomial 
Q is a factor of the product AB, then Q either divides A or Q divides B or there is 
a unique quaternion g G El and two quaternion polynomials Ai and Bi such that 
Q = {t — q){t — q), A = Ai{t — q) and B = {t — q)Bi. 

Proof. It is sufficient to prove the third case under the assumption gcd((5, Aq) = 
gcd((5, Bq) = I where Aq := mrpf A, Bq := mrpf B. There exist uniquely determined 
quaternions qi,... ,qr, gj,..., such that 

• Q = {t-qi){t-qj) = {t-ql){t-ql) for i G {!,.. .,r} and j G {!,... ,s}, 

• A = AqQa, B = BoQb where Qa ■= (t - qQ ■ ■ ■ (t - q 2 ){t - qi), Qb ■= 
{t - q[ )^- ^ 2 ) • • • (t - q'Qj_ and 

• gcd(AoAo, Q) = gcd(Po.Bo, Q) = 1. 

This follows from the quaternion version of the factorisation theorem for motion 
polynomials [5, Theorem 1]. 

As Q divides AB, we get that Q divides AqAABBq = AqAqQaQbBqBq. But 
then Q also divides QaQb- This can happen only if two neighbouring linear factors 
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are conjugate and, because of gcd((5, mrpf A) = gcd((5, mrpf B) = 1, implies qi = gj 
whence the claim follows with q ■= qi = D 


Proof of Theorem 1. Let T := 2xo + £(xii + X 2 j + a^sk) be the translation which 
moves a point x = xq + xA + a; 2 j + xsk to the origin 1 and let C = CT~^. Then 
the orbit of x with respect to C is equal to the orbit of the origin 1 with respect to 
C. Because we have deg(C') = deg(C') and the spherical degree defect of C equals 
the spherical degree defect of C, it is sufficient to prove the statement for C and the 
origin 1 (instead of C and x). From now on, we use the symbol C to abbreviate C. 
Let C = P + eQ, set G := mrpf P, m := deg G, P' := P/G and use (3) to compute 
a parametric equation y for the trajectory of the origin 1: 

y = PP+ 2PQ = G^P'P' + 2GP'Q = GP'P' + 2P'Q, 


where we write “=” for equality in the projective sense, modulo scalar (or real 
polynomial) multiplication. The trajectory’s degree is not larger than 2n — m. 

With N := gcd(G'P'P', mrpf(P'(5)) and r := deg the degree of the trajectory 
yisd=2n — m — r because we have y = {GP'P' + 2P'Q) /N. If P'P' /N and QQ/N 
are polynomials, the circularity is half the degree of 



GP'P' 

N 


P'P'QQ\ 

) 



QQ \ 

N )' 


This would then imply that the circularity is not less than n — m — r/2, which equals 
{d — to)/ 2 for d = 2n — m — r. Thus, we have to show that N divides both, P'P' 
and QQ. 

Moreover, it is sufficient to prove the claim only for a reduced Go ■= Pq + sQo as 
the claim in the non-reduced case follows easily. Set Nq := gcd(G'oPoPo, mrpf PqQo) 
where Gq := mrpf Pq and Pq := Pq/Gq. It is a useful fact that Nq has no real linear 
factor because such a factor would also divide Go (because Pg e has no real 
factor and hence PgPg has no real linear factor either) and Qo (because it divides 
PgQo and Pq has no real factor). But this contradicts the reducedness of Gq. Now 
we proceed by induction on the degree s of Nq. Note that s is zero or an even 
positive integer. 

In case of s = 0, Aq is a real constant and the claim is clear. 

Let s > 2. Then there is a monic quadratic irreducible real polynomial Mi which 
divides Nq. Set Ni := Nq/Mi. We claim that there always exists a quaternion q and 
two quaternion polynomials P{ and Qi such that Mi = {t — q){t — q), Pq = P[{t — q) 
and Qo = Qi{t — q). We handle the proof of this claim by distinguishing three cases; 

(1) Ml divides Gq. Then Mi can not divide mrpf Qo because of the reducedness 

of G. Since Mi is a divisor of iVo which divides PqQo, Mi must divide PqQo- 
Since Pq has no real polynomial factor, by Lemma I, there is a unique 
quaternion q and two polynomials P( and Qi such that Mi = {t — q){t — q), 
Po = P{{t - q) and Qo = Qi{t - q). _ 

(2) Ml does not divide Gq but divides mrpf Qo- Then Mi must divide PqPq. 
Since Pq has no real polynomial factor, by Lemma I, there is a unique 
quaternion q such that Mi = {t — q){t — q), Pq = P( {t — q). Then we set 
Qi '■= Qo{t — q)lMi which is a polynomial and satisfies Qo = Qi(^ — d)- 

(3) Ml divides neither Go nor mrpf Qq. By Lemma 1, there is a unique quater¬ 
nion q and two polynomials P{ and Qi such that Mi = {t — q){t — q), 
Pq = Piit - d) and Qo = Qi{t - q). 
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In all three cases, we have a new reduced motion polynomial Ci := Go-P( + eQi, 
such that 

Pq = Piit - q), Qo = Qi(t-q), gcd(GoPi'Pi',mrpf(P(Qi)) = ^o/Mi =: Ni. 

By induction hypothesis, Ni divides P(Pl and QiQi- By the derivation of P( and 
Qi, namely, Pq = P{{t — q) and Qo = Qi{t — q), we have that Nq divides PqPq and 
QoQo- This concludes the induction proof and also the proof of the theorem. □ 

Remark 2. Calling a rational motion generic if the primal part has no real factors 
one interpretation of Theorem 1 is as follows: A trajectory of a generic rational 
motion can be entirely circular (actually, only few exceptions are known). A similar 
statement for algebraic planar motions can be found in [2, Ch. XI]. There, also 
a more detailed discussion on the circularity of trajectories can be found. For 
planar rational curves, the geometric circularity conditions there imply the algebraic 
circularity conditions of Theorem 1. 

As one consequence of Theorem 1, among all rational curves of fixed degree, 
curves of high circularity can be generated by rational motions of low degree. In 
particular, we obtain the desired bound on the degree of a rational motion with a 
given rational trajectory: 

Corollary 1. If a rational curve of circularity c and degree d is a trajectory of the 
rational motion C = P + sQ, the degree of C is not less than d — c. If it is of degree 
d — c, the degree defect of the spherical motion component equals d — 2c. 

We will see below in Section 3.2 that the bound of Corollary 1 is sharp. 

Example 1. Let us illustrate above results by simple examples from literature. In case 
of degG = 1, the only generic rational motions are rotations with fixed axis. Their 
generic trajectories (trajectories that maximise the degree among all trajectories) 
are circles and, of course, entirely circular. Rational motions with degG = 2 are 
generated by reflecting a moving frame in one family of rulings on a quadric H, 
see [4]. The generic case is obtained if iL is a hyperboloid. In this case, generic 
trajectories are of degree four. Their entirely circularity has already been observed 
in [2, Ch. IX, §7]. The non-generic trajectories are circles. Points with non-generic 
trajectories lie on two skew lines that coincide if is a hyperboloid of revolution. 
Finally, if is a hyperbolic paraboloid, generic trajectories are just of degree three 
and circularity one [4]. 

3.2. Motion synthesis. Now we turn to the task of computing a rational motion 
C = P + eQ of minimal degree deg C = n with a given rational parameterised curve 
X = xq + xii + X 2 j + x^'k as trajectory. Denote the degree of the trajectory by d 
and its circularity by c. By Corollary 1 we have n > d — c. We will prove that, up 
to coordinate changes, exactly one solution of minimal degree n = d — c exists and 
we also provide a procedure for its computation. 

Theorem 2. Let x = a;o-l-a;ii-|-a: 2 j + a: 3 k be a reduced rational parametric equation 
of a rational curve such that x{oo) = 1. Then there exists a unique monic rational 
motion polynomial C = P + eQ of minimal degree, such that G(oo) = 1 and the 
trajectory of 1 is parameterised by x. 

A key ingredient in our proof of Theorem 2 is a version of the Euclidean algorithm 
to compute the left gcd of two polynomials F,G € This is a well-known concept 
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in the theory of polynomials over rings, see [10]. Call the polynomial L G IHI[t] a 
left factor or left divisor of F if there exists Q G IHI[t] such that F = LQ. The 
polynomial D G IHI[<:] is called left gcd of F and G, if U is a left divisor of F and G 
and any left divisor E of F and G also left divides D. The left gcd is unique up to 
right multiplication with a non-zero quaternion. 

The Euclidean algorithm in this context is based on polynomial right division. 
Given Ro,Ri G IHI[t], there exist Q 21 R 2 G IHI[1] such that Rq = R 1 Q 2 + R 2 and 
degi ?2 < degi?i. Note that the order of factors in the product R 1 Q 2 is important. 
If £) is a left divisor of i?o and i?i, then D is also a left divisor of R 1 Q 2 and of i? 2 - 
Conversely, a left divisor of R 2 and i?i also left divides Rq- Hence the left gcd of Rq 
and i?i equals the left gcd of i?i and i? 2 - Assuming degi?o > degi?i, we also have 
deg Rq > deg R 2 and we may recursively define sequences R 2 , R 3 , ■ ■. and Q 2 ,Q 3 , ■ ■ ■ 
of polynomials with strictly decreasing degree by 

Rk -2 = Rk-lQk + Rk, k>2. 

The recursion ends as soon as Rk = 0 whence the left gcd is D = Rk-i- Us¬ 
ing polynomial long division, an algorithmic implementation of this algorithm is 
straightforward. 

Lemma 2. If G G IHI[t] and M G ffi[f] is an irreducible (over Kj quadratic factor of 
GG, then M and G have a left gcd of positive degree. 

Proof. Use polynomial right division to compute Q,R G IHI[<] with G = MQ + R = 
QM + R and deg i? < 1. Then 

GC = {MQ + R){MQ + R) = M{MQQ + QR + RQ) + RR. 

If i? = 0, then M itself is a left gcd of G and M. Otherwise, M is a left factor of 
RR which implies M = RR. Hence, i? is a left factor of M and also of G. □ 

Lemma 3. For any polynomial G G IHI[t] of positive degree and a monic divisor 
R G M[t] of GG with gcd(i?, mrpf G) = 1, there exists a unique monic quaternion 
polynomial P (the left gcd of G and R) and a unique quaternion polynomial Q such 
that PQ = G and PP = R. 

Proof. The polynomial R has no real linear factor because such a factor necessarily is 
a divisor of C which contradicts gcd(i?, mrpf C) = 1. Using the Euclidean algorithm, 
we can compute the unique monic left gcd Pq of C and R, i.e., 

(4) PoQ = C, PoPi = R 

with Q, Pi G IHI[t]. We claim that Pq and Pi both have no real polynomial factor: 

• A real polynomial factor of Pq is also a real polynomial factor of G and R 
and contradicts gcd(i?, mrpf G) = 1. 

• A real polynomial factor of Pi must have a real quadratic factor M which 
is irreducible over M (because R has no real linear factor). That is, for some 
P[ G IHI[t] we have Pi = P[M. We can write R = PdPi = PoP{M = R'M 
where R' is real. From this we infer PqA' = PoPoP[. Because Pq has no 
real polynomial factor, the real polynomial PqPq divides R' and PqPqM 
divides R but also CG = QQPqPq (a multiple of R). This implies that M 
divides QQ and, by Lemma 2, Q and M have a left gcd P 2 with deg P 2 > 0. 
But then, by (4), P 0 P 2 is a left divisor of G and R and degPo-Pz > degP 2 . 
This contradicts the fact that Pq is left gcd of G and R. 
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Now we claim that Pi = Pq- Left multiplying the second equation in (4) with Pq 
we obtain PqR = PqPqPi. We have a real polynomial factor R on the left and 
a real polynomial factor PqPo on the right. Because neither Pq nor Pi have real 
polynomial factors, we have R = PqPq and Pi = Pq follows. Thus, with P := Pq, 
we have existence. 

As to uniqueness, we observe that P is necessarily the unique monic left gcd of 
C and R. Uniqueness of Q follows from uniqueness of P. □ 


Proof of Theorem 2. Denote the degree of the rational parametric equation x by 
d and its circularity by c. The condition a:(oo) = 1 implies degxo > degx^ for 
i = 1,2,3 and it is no loss of generality to assume that xq is monic. There exists a 
monic polynomial g of degree 2c and two relatively prime polynomials w (monic), 
y G M[t] with deg w = d — 2 c, deg y < 2 {d — c), such that 

(5) Xq = gw, x‘l + xl + xl = gy. 

Denote by C := P + eQ a (yet unknown) motion polynomial of minimal degree that 
parameterises the sought rational motion. The parametric trajectory of 1 equals x if 

(6) Xq + xii + X 2 S + xfk = PP + 2PQ where PQ + QP = 0. 

With D := a;ii + a; 2 j + a^ak, we have 

(7) DD = x\ + x\ + x\ = gy and wDD = wgy = XQy. 

By Lemma 3, there exist polynomials Pq G and Qo G such that 

^b-fo = 9 , PoQo = D. 


With P := wPq and Q := Qq/2 we then have 

PP = w^PqPq = w^g = wxq, 2PQ = wPqQq = wD = w{xii + X 2 j + xak). 

This shows that P and Q solve the first of the two equations in (6). The second 
equation follows from the vanishing of the scalar part of D: 0 = D+D = wD+wD = 
2{PQ + QP). Thus, the polynomial C = P + eQ with P, Q computed as above is 
one solution to our problem and it is of minimal degree n '■= d — c. 

Conversely, assume that C' = P' + eQ' is another solution of minimal degree 
deg C = n. Again, we have deg P' > deg Q' and P' is monic. By the first of the 
two equations in (6) we get 

(8) P'P' = w'xq, 2P'Q'= w'D = w'{xii +X 2 j +xqU), 


for some monic real polynomial w' of degree d — 2c. By (5), (7), and (8) we get 


W'^ W'^ W 

If we divide both sides by g and multiply with w', we get 

yw' = AwQ'Q'. 


As w and y are relatively prime, we get that w divides w'. But degw = degw' and 
w and w' are monic. Thus, we have 

(9) w = w' and y = AQ'Q'. 

Left-multiplying the second of two equations in (8) by Q' and using (9), we get 

, DQ' 2DQ' 

P = w -= = w -. 


2Q'Q' 


y 
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As w and y are relatively prime and w is monic, Pq '■= 2DQ'/y must be a monic 
polynomial. We have 


P' P' — 
^ 0 -* 0 ~ 


AQ'Q'DD DD 


yZ 


= - =g, 2P'Q' = 

y 


^Q'Q'D 

y 


= D. 


Because of degPg = degPg and because the left gcd of g and D is unique, this 
implies Pg = Pg and Q' = Qo/2 = Q. This concludes the proof of uniqueness. □ 


Remark 3. In case of a rational curve of circularity zero, the primal part P of the 
rational motion C = P + eQ is a real polynomial and the trivial solution consists 
of the curvilinear translation along the given curve. In all other cases. Theorem 2 
guarantees solutions of lesser degree than curvilinear translations. 


Our formulation of Theorem 2 requires a special moving point, the origin 1, to 
generate the rational curve and produces a motion polynomial C with the special 
property 0 ( 00 ) = 1. These conditions are just coordinate dependencies and are 
convenient for our algebraic formulation and proof. If we allow coordinate changes, 
we may translate 1 to an arbitrary moving point (this amounts to replacing C by 
CT~^ with a suitable translation, as in the proof of Theorem 1) and orient the 
coordinate axis arbitrarily (this amounts to right-multiplying C with r G H). Thus, 
we can re-formulate Theorem 2 as 


Corollary 2. There is a unique (up to coordinate changes) rational motion of 
minimal degree in the dual quaternion model of rigid body displacements with a 
prescribed rational trajectory. If the trajectory is of degree d and circularity c, the 
minimal motion is of degree n = d — c. 

Remark 4. Frequently, the degree of a rational motion is defined as maximal degree 
of one of its trajectories. With this concept of degree. Corollary 2 is no longer true. 
A counter example is a circle which occurs as trajectory of a circular translation and 
of a rotation about its centre. For both motions the maximal degree of a trajectory 
is two. 


Let us illustrate Theorem 2 by two examples: 

Example 2. The rational curve given by the parameterised equations 

Xo = + 2t + 2){t^ + 2t + 5)(^ + 1), Xi = 2{t^ - 5)(t^ + 2i + 2), 

X2 = —4(t + 5)(t + 1)^, X 3 = — 2t{t + 2)(t + 5)(i + 1) 
is of degree five. From the factorisation 

x\ + x\ + x\ = 8 (t^ -I- 2t -I- 5)(t^ -I- 4t -I- 5)(t^ -I- 2t -I- 2)^ 

we see that its circularity is two. By Theorem 2, it is trajectory of a rational motion 
of minimal degree three whose primal part, by Theorem 1, has precisely one linear 
factor. With y, w and g defined by 

g = (t^ 2t 2)(t^ 4 “ 2t 4-5), 'u: = t4“l, y ~ 8 (t^ 4- 4t 4- 5)(t^ 4 - 2t 4 - 2), 

we have 

Xo = gw, xl + xl + xl = gy. 

Setting D = xii 4- X2i 4- xak and using the Euclidean algorithm we compute, 
according to Lemma 3, the left gcd Pg of D and y. With 

Po = t^ + 2 t+l-{t+l)i+ {2t + 2)j - 2k 
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and 

C^Q = —2t — 2 — (2i^ -t- At 2)i -t- (2t -t- 6)j -t- (2t^ -t- St -t- 6)k 
we have g = PqPo, D = PqQq. Then we set P = wPq and Q = Qo/2. As expected, 
P has a real polynomial factor of degree 1. The motion polynomial we want is 
C:=P + eQ. 

Example 3. The Viviani curve [11] is given by 

xo = (l + t^)^, xi = (1 - t)^(l + X2 = 2 t{l - t){l + t), xs = 2 t{l + t"^). 

In order to meet the requirements of Theorem 2, we translate it by the vector 
(—1,0,0) and obtain 

(10) xo = (l + t^)^, xi =-4t^, X2 = 2t(l - t)(l + t), xs = 2 t{l + t '^)■ 

This curve lies on the unit sphere and is entirely circular. As in the previous example, 
we compute 

g = gcd{xo,xl + xl + xl) = {1 + t'^)'^, w = xo/g=l, y = {xj + xl + xD/g = St'^. 
The left gcd of D := xi'i + X 2 j + xsk and g is 

P^ = f- t{j + k) - i, 
the right quotient of D and Pq is 

Qo = 2t(k-j), 

and the minimal motion to the curve (10) is 

C = - (j + k - e(j + k))t - i. 

We can simplify by conjugating with T := 1 + f ei whence we get 

C' := TCP = f- t{i + k)-i={t- k)(t - j). 

This shows that C is the composition of two rotations about the second coordinate 
axis and the third coordinate axis with equal angular velocities, i.e., the motion 
generated by the rolling of a spherical circle of radius ir/A on a spherical circle of 
the same radius. This is illustrated in Figure 1. There, the moving frame is rigidly 
attached to the rolling circle. 

4. Discussion of results 

This article unveiled some relations between rational motions and their trajectories. 
A rational curve occurs as trajectory of a unique (up to coordinate changes) rational 
motion of minimal degree in the dual quaternion model. This was a surprise to the 
authors as a mere trajectory seems to leave a lot of freedom for the construction of 
a suitable rational motion. Apparently, the requirement for minimality is rather 
restrictive. 

Calling a rational curve “generic” if its circularity is zero and a rational motion 
“generic” if its primal part has no real factors (the spherical motion component has 
full degree), we may also say that rational motions of minimal quaternion degree 
to a generic trajectory are non-generic. Conversely, the trajectories of a generic 
rational motion are entirely circular and hence non generic. 

In conjunction with the factorisation of generic rational motions [5] and its 
extension [8] to non-generic motion polynomials, our results contribute to a variant 
of Kempe’s Universality Theorem [3, Section 3.2] for rational space curves. Via 
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Figure 1. Minimal motion of the Viviani curve (rolling of spherical circles) 


motion factorisation, it is possible to construct linkages to generate a given rational 
motion and hence also a given rational trajectory. The dual quaternion degree of 
the generating motion is directly related to the number of links and joints in the 
mechanism. As a consequence, rational curves can be generated by spatial linkages 
with much fewer links and joints than those implied by the asymptotic bounds for 
algebraic space curves given by [1]. For circular curves, the numbers of links and 
joints are even less. A precise formulation and a rigorous proof will be worked out 
in a future paper. 

Our main result in Corollary 2 raises questions. Is a similar uniqueness statement 
true for algebraic curves and algebraic motions? What rational motions are not 
minimal for any of their trajectories? The Darboux motion and the circular transla¬ 
tion of Remark 1 are examples. Finally, determine simple (classes of) curves with 
simple minimal motion as in Example 3 and exploit their low degree and rationality 
in a CAGD or engineering context. 
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